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The  paper  aims  to  overview  various  mathematical  models  for  modeling  the  simultaneous  heat  and  mass 
transfer  process  in  the  liquid  desiccant  dehumidifier.  Firstly,  the  dehumidification  principle  is  introduced 
briefly.  Then  the  models  are  interpreted  in  terms  of  two  classes  of  dehumidifiers.  For  the  adiabatic 
dehumidifier,  the  models  are  mainly  classified  into  three  types:  finite  difference  model,  effectiveness  NTU 
( e-NTU )  model,  and  simplified  models.  For  the  internally  cooled  dehumidifier,  there  are  also  three  kinds 
of  models:  models  without  considering  liquid  film  thickness,  models  considering  uniform  liquid  film 
thickness,  and  models  considering  variable  liquid  film  thickness.  This  review  is  meaningful  for 
comprehending  the  development  process  and  research  status  of  the  models  and  choosing  suitable 
models  for  prediction.  In  addition,  some  suggestions  are  proposed  for  the  model  improvement. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


Contents 


1.  Introduction . 587 

2.  Problem  description . 588 

2.1.  The  mechanism  of  heat  and  mass  transfer  in  liquid  desiccant  dehumidifier . 588 

2.2.  Structure  of  the  dehumidifier . 589 

3.  Models  for  adiabatic  dehumidifier . 589 

3.1.  Finite  difference  model . 589 

3.2.  Effectiveness  NTU  {e-NTU)  model . 591 

3.3.  The  simplified  models . 592 

3.4.  Summary . 593 

4.  Models  for  internally  cooled  dehumidifier . 594 

4.1.  Models  without  considering  liquid  film  thickness . 594 

4.2.  Models  considering  uniform  liquid  film  thickness . 595 

4.3.  Models  considering  variable  liquid  film  thickness . 596 

4.4.  Summary . 597 

5.  Conclusion . 597 

Acknowledgments . 598 

References . 598 


1.  Introduction 
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E-mail  address:  hong-xing.yang@polyu.edu.hk  (H.  Yang). 


With  the  acceleration  of  urbanization  and  improvement  of 
people'  living  standard,  a  larger  proportion  of  building  energy 
consumption  will  be  needed  to  keep  a  comfortable  indoor  envir¬ 
onment  [1  ].  But  it  is  well-known  that  the  traditional  air  condition¬ 
ing  system  is  notorious  as  a  result  of  heavy  dependence  on  electric 
power,  limited  ability  of  humidity  control,  and  occurrence  of  wet 
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Nomenclature 

A  specific  surface  area  per  unit  volume  [m-1] 

CP  specific  heat  []  kg-1  K-1] 

Csat  saturation  specific  heat  [J  kg-1  I<-1] 

dh  hydraulic  diameter  [m] 

Dm  mass  diffusion  coefficient  [m2  s-1] 

Dt  thermal  diffusivity  [m2  s-1] 

g  gravity  [ms-2] 

G  specific  mass  flow  rate  [kg  m-2  s-1] 

G'  mass  flow  rate  [kg  s-1] 

h  specific  enthalpy  [kj/kg] 

he  enthalpy  of  humid  air  in  equilibrium  with  liquid 

desiccant  [kj  kg] 

he,eff  effective  saturation  enthalpy  [kj  kg],  Eq.  (18) 

H  height  of  the  dehumidifier  [mj 

k  heat  conduction  coefficient  [W  m-1  K-1] 

L  length  of  the  dehumidifier  [m] 

Le  Lewis  number 

m  water  condensation  rate  [g  s  1  ]  or  per  unit  cross- 

sectional  area  [g  m-2  s-1],  Eq.  (23) 
m*  capacity  ratio  similar  to  the  one  used  in  sensible  heat 
exchangers 

JVf  molecular  weight  [g  mole-1] 

NTU  number  of  transfer  units 

Nu  Nusselt  number  (dimensionless) 

P  pressure  [Pa] 

Pa  partial  vapor  pressure  in  air  [Pa] 

Ps  partial  vapor  pressure  over  the  solution  [Pa] 

Pt  total  pressure  [Pa] 

P  dimensionless  vapor  pressure  difference  ratio 

Q  heat  transferred  from  solution  to  water  [kW  m_1] 

T  temperature  [K] 

T  dimensionless  temperature  difference  ratio 

u  velocity  [m  s- ’] 

V  volume  [m3] 

w  humidity  ratio  [kg  H20  kg-1  dry  air]  or  the  width  the 

dehumidifier  [m] 


We  humidity  ratio  of  humid  air  in  equilibrium  with  liquid 
desiccant  [kmol  H20  (kmol  air)-1] 

Weeff  effective  humidity  ratio  [kmol  H20  (kmol  air)-1], 

Eq.  (19) 

X  desiccant  concentration  [kg  desiccant  kg-1solution] 

Xw  concentration  of  water  in  solution 
[kg  water  kg-1solution] 

Xv  concentration  of  water  vapor  in  air  [kgH20  kg-1] 

Greek  letters 

ac  heat  transfer  coefficient  [W  m-2  K-1] 

a'c  heat  transfer  coefficient  corrected  for  simultaneous 

mass  and  heat  transfer  [W  m-2  K-1] 
aD  mass  transfer  coefficient  [kg  m-2  S-1] 
a’D  molar  mass  transfer  coefficient  [kmol  m-2  S-1] 

S  dimensionless  temperatures  ( T-Tr)lh ,  Eq.  (33) 

p  density  [kg  m-3] 

p  dynamic  viscosity  [  kg  m  - 1  s  - 1  ] 

v  kinetic  viscosity  [m2  s-1] 

e  air  side  effectiveness 

eHE  heat  exchanger  effectiveness 

X  latent  heat  of  vaporization  [kj  kg-1] 

5  film  thickness  [m] 

A  change  of  or  difference  between  parameters 

Subscripts 

a  air 

c  critical 

/  cooling  fluid,  like  water,  air,  refrigerant 

i  inlet 

inf  interface 

o  outlet 

p  primary  air 

r  secondary  (return  air) 

s  desiccant  solution 

v  water  vapor 


surface  for  breeding  mildew  and  bacteria  and  so  on  [2],  Thus,  to 
reduce  the  energy  consumption  in  buildings  and  improve  the 
indoor  air  quality,  the  liquid  desiccant  assisted  air  conditioning 
system  has  drawn  more  and  more  attention  [3-7], 

The  major  component  of  interest  regarding  heat  and  mass 
transfer  of  such  a  system  is  the  dehumidifier.  Compared  with  the 
experimental  research,  the  simulation  method  is  more  time  and 
cost  saving.  Also,  some  parameters  in  the  dehumidifier  interior  can 
be  observed  by  simulation  while  it  is  impossible  to  be  achieved  by 
experiment.  Most  importantly,  the  verified  simulation  models  are 
effective  tools  to  assess  and  optimize  similar  dehumidifiers.  There¬ 
fore,  a  large  amount  of  studies  have  been  done  to  establish 
reasonable  mathematical  models  for  evaluating  the  liquid  desic¬ 
cant  dehumidifiers.  However,  there  is  short  of  detailed  and  specific 
summary  of  the  models  until  now.  Thus,  it  is  meaningful  to  classify 
and  assess  the  models,  which  will  provide  useful  suggestion  for 
future  research. 

In  the  paper,  the  function  principle  of  the  liquid  desiccant 
dehumidifier  is  introduced  firstly.  Based  on  whether  there  is  heat 
removal,  the  dehumidifier  is  divided  into  two  types:  adiabatic  and 
internally  cooled  dehumidifier.  Correspondingly,  the  models  are 
summarized  in  two  respects  in  terms  of  the  different  structures. 
For  each  model,  the  assumptions,  governing  equations,  boundary 


conditions  and  other  relevant  information  are  provided.  The 
applied  conditions,  development  process,  and  research  status  of 
the  simulation  models  are  also  presented.  In  addition,  some 
suggestions  are  put  forward  for  the  model  improvement. 


2.  Problem  description 

2.1.  The  mechanism  of  heat  and  mass  transfer  in  liquid  desiccant 
dehumidifier 

It  is  well  know  that  in  the  dehumidifier,  complicated  heat  and 
mass  transfer  occurs.  The  driving  force  for  heat  transfer  is  the 
temperature  difference  between  the  air  and  desiccant  solution, 
and  for  mass  transfer  is  the  water  vapor  pressure  difference 
between  the  air  and  the  surface  of  the  desiccant  solution.  The 
most  classic  and  typical  mass  transfer  theories  include  film  theory, 
penetration  theory  and  surface  renewal  theory.  The  theory  used 
most  for  the  dehumidifier  is  the  film  theory.  It  is  Nernset  [8]  who 
proposed  the  film  theory  first  in  1904.  He  assumed  that  the  whole 
resistance  of  mass  transfer  in  a  given  phase  lied  in  a  thin  and 
stagnant  region  of  that  phase  at  the  interface.  This  region  is  called 
film.  Based  on  it,  Whiteman  [9]  developed  the  two-film  theory. 
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The  specific  transport  mechanism  is  shown  in  Fig.  1.  Where,  PB  is 
the  partial  pressure  of  component  B  in  the  gas  phase  and  XB  is  the 
mole  fraction  of  component  B  in  the  liquid  phase. 

The  two  film  theory  is  very  easy  to  understand  and  apply,  but  it 
has  several  drawbacks.  Firstly,  it  is  not  reasonable  to  predict  that 
the  rate  of  mass  transfer  is  directly  proportional  to  the  molecular 
diffusivity.  Secondly,  it  is  difficult  to  decide  the  thickness  of  the 
two  laminar  sub  layers  by  experiment.  Finally,  the  convective  mass 
transfer  in  the  thin  films  is  neglected,  so  the  theory  is  only  suitable 
for  the  steady  mass  transfer  process. 

In  the  dehumidification  process,  some  quantity  of  heat  is  given 
out  as  well,  including  the  phase  change  heat  and  dilution  heat.  In 
all  of  the  dehumidifier  models,  the  dilution  heat  is  neglected  as  it 
is  much  smaller  compared  with  the  phase  change  heat  of  water 
vapor  [10], 

2.2.  Structure  of  the  dehumidifier 

As  for  the  structure  of  the  dehumidifier,  according  to  whether 
there  is  heat  output,  the  dehumidifiers  can  be  classified  into 
adiabatic  and  internally  cooled  dehumidifier  [11].  The  diagrams 
of  two  dehumidifier  structures  are  shown  in  Fig.  2. 

In  an  adiabatic  dehumidifier,  the  air  and  liquid  desiccant 
contact  directly  with  each  other.  In  the  early  stage,  the  research 
was  concentrated  on  the  structure  of  spray  tower  as  a  result  of  its 
simple  construction  and  large  specific  surface  area  [12],  However, 
in  the  spray  tower,  the  desiccant  solution  is  generally  broken  into 
small  droplets,  so  the  problem  is  sometimes  serious  of  mist 
generation  and  carryover  of  liquid  droplets  in  the  air  stream.  Then, 
the  packed  tower  was  used  widely  because  it  is  more  compact  and 
can  provide  a  higher  residence  time,  lower  liquid  pressure  loss  and 


Interface 


Mass  Transfer  Direction 
Fig.  1.  Schematic  diagram  of  two-film  theory. 


lessen  the  carryover  problem.  In  1980,  Factor  and  Grossman 
verified  the  possibility  of  employing  the  packed  tower  as  dehu¬ 
midifier  by  theoretical  analysis  and  experiment  [13].  As  for  the 
padding  materials,  the  random  packing  like  ceramic  [14],  plastic 
and  polypropylene  pall  ring  are  popular  first  [15,16],  Then  some 
structured  packing  materials  are  employed  to  optimize  the  flow 
and  reduce  the  resistance  in  the  dehumidifier,  such  as  the  stainless 
steel  corrugated  orifice  plate  [17],  celdek  [18,19]  and  so  on. 

In  the  adiabatic  dehumidifier,  the  temperature  rise  of  the 
desiccant,  resulted  from  the  latent  heat,  worsens  its  dehumidifica¬ 
tion  performance,  thus  its  dehumidification  efficiency  is  relatively 
lower.  A  solution  is  to  increase  the  desiccant  flow  rates  to  achieve 
good  dehumidification  levels.  However,  the  high  desiccant  flow 
rates  and  the  followed  higher  flow  rates  of  the  regenerated 
desiccant  solution  reduce  the  coefficient  of  performance  of  the 
liquid  desiccant  cycle  [20,21  ].  In  addition,  the  desiccant  particles 
are  much  easier  to  be  entrained  by  the  air  and  therefore  pollute 
the  indoor  environment. 

To  solve  the  above  problems,  the  internally  cooled  dehumidifier 
was  developed.  In  an  internally  cooled  dehumidifier,  besides  the 
contact  between  air  and  desiccant,  some  cold  source  which  can 
provide  cool  fluid  like  air  or  water  is  added  to  take  away  the  latent 
heat  produced  in  the  process  of  dehumidification,  which  can  be 
regarded  as  an  isothermal  process  in  general.  The  internally  cooled 
dehumidifier  has  been  popular  since  the  1990s  [22,23],  As  the 
latent  heat  is  removed  from  the  dehumidifier,  it  reduces  the 
temperature  rise  of  the  solution  and  air,  resulting  in  efficiency 
improvement  [24].  Meanwhile,  it  allows  lower  desiccant  flow  rates 
in  the  internally  cooled  dehumidifier  so  as  to  reduce  the  pollution 
problem.  But  the  internally  cooled  dehumidifier  has  more  com¬ 
plicated  structure  than  the  adiabatic  one.  For  example,  to  increase 
the  contact  area,  a  fin  structure  is  widely  used  in  the  internally 
cooled  dehumidifier  or  other  heat  and  mass  transfer  devices 
[25-29], 

3.  Models  for  adiabatic  dehumidifier 

There  are  mainly  three  types  of  mathematical  models,  includ¬ 
ing  the  finite  difference  model,  effectiveness  NTU  (e-NTU)  model 
and  the  simplified  solutions. 

3.1.  Finite  difference  model 

In  1980,  Factor  et  al.  [13]  promoted  a  theoretical  model  to 
predict  the  performance  of  a  countercurrent  packed  column  air- 
liquid  contractor,  based  on  the  model  for  adiabatic  gas  absorption 
put  forward  by  Treybal  in  1969.  In  the  model,  the  whole  dehumi¬ 
difier  is  divided  into  n  parts,  as  shown  in  Fig.  3. 

To  simplify  the  complexity  of  the  heat  and  mass  transfer 
process,  several  assumptions  were  made:  (1)  the  flow  of  air  and 
desiccant  were  assumed  as  the  slug  flow,  (2)  the  process  was 
adiabatic,  (3)  the  properties  of  the  gas  and  liquid  were  assumed 


Fig.  2.  The  structure  diagram  of  two  dehumidifiers. 
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constant  across  the  differential  element,  which  meant  the  gradi¬ 
ents  only  exist  at  the  z  direction,  (4)  both  of  the  heat  and  moisture 
transfer  areas  were  equal  to  the  specific  surface  area  of  the 
packing,  (5)  it  was  negligible  of  the  non-uniformity  of  the  air 
and  solution  flows,  (6)  in  the  flow  direction,  no  heat  and  moisture 
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Fig.  3.  Heat  and  moisture  exchange  model  in  the  countercurrent  adiabatic 
dehumidifier. 


transfer  occurred,  (7)  the  resistance  to  heat  transfer  in  the  liquid 
phase  was  negligible,  and  (8)  the  interface  temperature  was  equal 
to  the  bulk  liquid  temperature.  Based  on  the  above  assumptions, 
the  main  governing  equations  were  stated. 

According  to  the  mass  balance  in  the  control  volume, 

dGs  =  GadW  (1) 


According  to  the  interface  mass  and  sensible  heat  transfer  rates, 
the  air  humidity  change  was, 


dW  a'DMvA  (\-Ps/Pt\ 
dz  Ga  \l-Pa/PtJ 


(2) 


According  to  the  interface  sensible  heat  transfer  from  the  air  to 
solution  side  and  the  energy  balance  on  the  gas  side,  the  air 
temperature  change  was, 

dT a  _  a'C  a A(  Ta  —  Ts) 

~dz~  GaCp# 


/  &  _ GaCPy(dW/dz) _ 

c’a  1-exp  [GaCPy(dW/dz)/(ac,aA)] 

where  ac,a  and  a'Ca  are  the  heat  transfer  coefficient  (sensible)  of 
the  gas  side  and  that  coefficient  corrected  for  simultaneous  mass 
and  heat  transfer  by  applying  the  Ackermann  correction,  which  is 
one  method  to  take  into  account  the  effect  of  mass  transfer  on  the 
temperature  profile  with  an  Ackermann  correction  factor. 

Finally,  the  boundary  conditions  were:  z  =  0  Ts=Tsj,  Gs=Gsi, 
X=Xr,  z  =  H,  Ta=Tai,  Ga=Gai,  W=Wi. 

Since  the  above  differential  equations  cannot  be  solved  analy¬ 
tically,  the  most  basic  solution  is  a  numerical  integration  along  the 
height  of  the  dehumidifier.  To  begin  the  calculation,  one  end  of  the 
dehumidifier  must  be  chosen  as  the  start  point.  For  the  counter- 
current  flow  pattern,  it  needs  to  presume  the  outlet  conditions  for 
one  of  the  fluids.  Solving  the  above  equations  from  the  bottom  to 
the  top  of  the  dehumidifier,  with  the  boundary  conditions,  a  set  of 
calculated  inlet  solution  parameters  are  obtained.  By  comparing 
the  calculation  results  with  the  real  values,  the  supposed  existing 


Fig.  4.  Calculation  flow  chart  of  countercurrent  pattern. 
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Fig.  5.  Schematic  of  the  cross  flow  dehumidifier/regenerator  [39]. 
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Fig.  6.  A  two  dimensional  schematic  of  the  cross  flow  dehumidifier/regenerator 

[39]. 


solution  variables  are  adjusted.  The  calculation  will  last  until  the 
final  results  are  very  close  to  the  real  values.  And  the  general 
calculation  flowchart  for  the  finite  difference  model  of  the  counter 
flow  pattern  can  be  summarized  in  Fig.  4. 

In  the  later  study,  Gandhidasan  et  al.  [30]  utilized  the  similar 
model  to  study  various  parameters  on  the  packing  height  of 
the  packed  tower.  Then,  Lazzarin  et  al.  [31]  gave  more  specific 
explanation  of  the  calculation  method  in  Appendix  A  of  the  litera¬ 
ture.  Oberg  and  Goswami  [32]  applied  a  finite  difference  model 
similar  to  Factor  and  Grossman's  to  verify  the  experimental  results. 
Taking  account  the  insufficiently  wetted  packing  and  the  different 
factor  when  transfer  the  fc-type  mass  transfer  coefficients  to  the  F- 
type  one,  Fumo  and  Goswami  [33]  improved  Oberg  and  Goswami's 
mathematical  model  by  modifying  the  transfer  surface.  In  addition,  a 
correction  factor  CF  was  introduced  to  modify  the  correlation  of  the 
wetting  surface.  By  comparing  the  results  of  simulation  and  experi¬ 
ment,  it  was  found  that  the  calculation  results  of  the  adapted  model 
agreed  well  with  experimental  results. 

All  of  the  above  models  introduced  a  coefficient  a'Ca  to  describe 
the  simultaneous  heat  and  mass  transfer  with  the  Ackermann 
correction.  Khan  and  Ball  [34]  promoted  another  solution  to  deal 
with  the  simultaneous  transfer  process.  Both  heat  and  mass 
transfer  processes  were  assumed  to  be  gas  controlled,  so  the 
interface  temperature  was  the  temperature  of  the  bulk  liquid 
and  the  heat  transfer  rate  across  the  air  film  from  the  bulk  air  to 
the  interface  was  equal  to  that  entering  the  liquid  side, 


(5) 


Similarly,  the  mass  transfer  across  the  interface  was  equal  to 
the  change  in  humidity  ratio, 

CadW  =  aDaA(We  —  W)dz  (6) 

Then,  the  humid  air  specific  enthalpy  change  could  be  written 
as, 


dha  —  CpadTa  +  dW  ■  \CpV(T a  —  Tr)  +  2]  (7) 

By  substituting  Eqs.  (6)  and  (7)  to  (5),  the  air  enthalpy  change 
along  the  flow  direction  was  obtained.  Here  it  is  rewritten  in  a 
simpler  way  appeared  in  the  later  literature, 


dha 

dz 


NTU  ■  Le 
H 


(he  —  /ia)+2 


■  (We  -  W ) 


(8) 


In  the  above  equations,  Le  and  NTU  were  defined  as 


Le  = 


«c 

al)Upa 


(9) 


NTU  = 


aDAV 


(10) 


In  this  way,  the  coupled  heat  and  mass  transfer  were  consid¬ 
ered  together.  In  the  following  research,  this  handling  method  is 
more  popular  than  the  Ackermann  correction. 

The  finite  difference  model  has  been  widely  used  for  the 
countercurrent  dehumidifier  [35-38],  For  cross  flow  configuration, 
Liu  et  al.  [39]  proposed  a  model  for  the  heat  and  mass  transfer 
process  in  a  cross  flow  adiabatic  liquid  desiccant  dehumidifier/ 
regenerator.  The  physical  and  mathematical  models  are  described 
in  Figs.  5  and  6,  respectively. 

The  governing  equations  of  energy,  water  content,  and  solute 
mass  balances  in  a  differential  element  were, 


G'a  _  dha  1  d(G'hs)  _  q 
H  dz  L  dx  ~ 


G«  W  1  ^  =  0 

H  dz  L  dx 


(12) 


d(G'  X)  =  0 


(13) 


The  energy  and  mass  transfer  in  the  interface  of  the  air  and 
desiccant  solution  were  expressed  in  Eq.  (8)  and  the  following 
Eq.  (14), 


dW 

dz 


NTU 

L 


■( We-W ) 


(14) 


Like  some  other  papers,  Le  was  supposed  to  be  one  in  the 
model.  However,  the  value  of  NTU  was  correlated  based  on  the 
corresponding  experimental  data  in  the  paper. 

Niu  [40]  also  established  a  two  dimensional  mathematical 
model  for  the  cross-flow  adiabatic  dehumidifier,  and  the  mass 
transfer  coefficient  was  gained  from  the  experimental  data.  Woods 
and  Kozubal  [41]  applied  the  finite  difference  model  to  study  the 
performance  of  a  desiccant-enhanced  evaporative  air  conditioner 
and  the  simulation  results  showed  good  agreement  with  the 
experimental  ones. 


3.2.  Effectiveness  NTU  (e-NTU)  model 

Stevens  et  al.  [42]  reported  an  effective  model  for  liquid- 
desiccant  heat  and  mass  exchanger,  which  was  developed  from  a 
computationally  simple  effectiveness  model  for  cooling  towers 
[43],  Except  for  the  assumptions  of  the  finite  difference  model,  two 
additional  assumptions  were  included.  One  was  the  assumption  of 
the  linear  relationship  of  saturation  enthalpy  and  temperature,  the 
other  one  was  the  neglect  of  the  water  loss  term  for  the  solution 
energy  balance.  In  addition,  an  ‘effective’  heat  and  mass  transfer 
process  was  assumed. 


GaCp,adTa  —  ac^aA(Ts  —  Ta)dz 
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The  main  equations  and  calculation  process  of  e-NTU  model 
were  summarized  as  follows, 


(1)  Calculated  the  Number  of  Transfer  Units  (NTU)  by  Eq.  (10). 

(2)  In  terms  of  the  similarity  with  the  heat  exchanger,  the 
effectiveness  of  the  countercurrent  flow  dehumidifier  could 
be  expressed  as 

|  _  g-NTU(l-m*) 

£=  l  _m*e-NTU(l-m*) 


where  m*  was  a  capacitance  ratio,  defined  analogous  to  the 
capacitance  ratio  used  in  sensible  heat  exchangers,  and  it  was 
given  in  the  following  equations, 


GaCsat 

s 


(16) 


where  Csat  was  the  saturation  specific  heat,  and 
Csar  —  (dhe/dTs). 

(3)  With  NTU  and  e,  the  air  outlet  enthalpy  could  be  obtained  with 
the  following  equation, 


dao  —  haj  +  e(h(:  —  hai ) 


(17) 


(4)  Used  an  energy  balance  to  calculate  the  solution  outlet 
enthalpy. 

(5)  Then  an  ‘effective’  saturation  enthalpy  was  found  by 


he.eff  —  da  i  -f- 


(Tq.o  haj 

\-e-NTU 


(18) 


(6)  Using  the  enthalpy  and  saturated  conditions,  the  effective 
humidity  ratio  Ye,eff  could  be  obtained. 

(7)  Then,  by  the  following  equation,  the  air  outlet  humidity  ratio 
could  be  calculated, 

W0  =  We,eff+(Wi-W^eff)e-NTU  (19) 


(8)  With  the  mass  balance  and  known  inlet  and  outlet  parameters, 
all  of  the  outlet  parameters  were  acquired. 

In  the  later  study,  Sadasivam  and  Balakrishnan  [44]  pointed  out 
that  the  definition  of  NTU  based  on  the  gas  mass  velocity  in 
Stevens's  model  was  not  appropriate  when  the  minimum  flow 
capacity  was  the  liquid  [45],  Thus,  the  gas  mass  velocity  C'a  in 
Eq.  (10)  was  changed  to  the  minimal  mass  velocity  of  gas  and 
liquid. 

As  for  the  e-NTU  model,  there  are  much  fewer  literatures  than 
the  finite  difference  model.  In  the  following  study,  Ren  [46] 
developed  the  analytical  expressions  for  the  e-NTU  model  with 
perturbation  technique.  The  model  accounted  for  the  nonlinea¬ 
rities  of  air  humidity  ratio  and  enthalpy  in  equilibrium  with 
solutions,  the  water  loss  of  evaporation  and  the  variation  of  the 
solution  specific  heat  capacity. 

3.3.  The  simplified  models 

It  can  be  found  that  both  the  finite  difference  model  and  e-NTU 
model  require  numerical  and  iterative  computations.  Thus,  both  of 
them  are  not  suitable  for  hourly  performance  evaluation. 

Khan  and  Ball  [34]  developed  a  simplified  algebraic  model. 
With  the  finite  difference  model,  about  1700  groups  of  data  were 
analyzed.  The  following  functions  were  deduced, 

W„=  n0  +  n  i  W,- + n2  Ts ,- + n3  Tj  j  (20) 

Wo  =  m0  +  m-lWi+m2TatO+m3Ta0  (21) 


The  above  two  equations  can  be  employed  to  predict  the  exit 
air  temperature  and  humidity  ratio  easily.  However,  as  they  were 
fitted  based  on  the  data  with  certain  operating  solution  concen¬ 
tration  and  solution  to  air  mass  ratio,  they  might  not  be  suitable 
for  other  conditions. 

Liu  et  al.  [47]  fitted  some  empirical  correlations  to  estimate 
the  performance  of  a  cross-flow  or  counter-flow  liquid  desiccant 
dehumidifier.  The  essence  of  the  methodology  is  obtaining  the 
empirical  expression  of  enthalpy  and  moisture  effectiveness  by 
experiment. 

Gandhidasan  [48]  reported  a  very  simple  analytical  solution  to 
predicate  the  rate  of  moisture  removal  for  the  dehumidification 
process.  Referred  to  previous  work  [49],  the  author  promoted 
the  dimensionless  moisture  and  temperature  difference  ratios. 
By  combining  the  aforementioned  two  definitions,  the  energy 
balance  equation  could  be  expressed  as  follows, 


Cp.aT(TaJ  -  rs,i)  +  ^  ■  ^P(Pa.i  -  Ps.d  =  -p,(Ts.o  -  TsJ) 


(22) 


According  to  the  literature  [50],  the  relationship  between  the 
rate  of  moisture  removal  m  and  the  partial  pressure  of  water  vapor 
could  be  deduced  as 


Ps.i 


mPt  Mi j 
GaP-Mv 


(23) 


In  addition,  the  desiccant  outlet  temperature  was  easily  calcu¬ 
lated,  as  shown  in  the  equation, 


Ts.o  = 


Tsj  —  ^heTcj 
(1  -che) 


(24) 


Finally,  substituting  Eqs.  (23)  and  (24)  to  (22),  the  moisture 
removal  rate  m  was  given  as, 


1  G'sCpseHE 

2  .(1  -che) 


{Ps.i  —  Tc.i)  —  GaCp  aT(Tai  —  TSj) 


(25) 


The  method  was  rather  simple  yet  it  involved  lots  of  assump¬ 
tions  and  limitations.  Except  some  common  assumptions,  it  also 
required  that  the  desiccant  inlet  temperature  was  different  from 
the  air  inlet  temperature,  and  the  desiccant  temperature  leaving 
the  dehumidifier  was  equal  to  that  of  the  desiccant  entering  the 
heat  exchanger. 

Chen  et  al.  [51]  constructed  an  analytical  model  on  the  basis  of 
the  finite  difference  model  for  countercurrent  and  concurrent  flow 
pattern.  The  physical  model  is  similar  to  that  of  Fig.  3.  Firstly,  a 
mathematical  model  was  built  following  the  model  promoted  by 
Khan  and  Ball  [34], 

Then,  two  parameters  were  introduced  for  derivation  conve¬ 
nience, 


K a  —  Le  ■  C p  a  l  a  +  7.  •  W 


(26) 


Ke  —  he  ■  Cp a  ■  Ts  +2  ■  IV c 


(27) 


By  combining  the  mass,  energy  conservation  equations,  mass 
and  energy  transfer  equations  at  the  interface,  the  change  of  Ke 
along  the  flow  direction  was, 


dKe  „NTU 

~r1=m —n- 

dz  H 


lja  i+rrf^Ke  ~  +  (Ec  —  1  )Cp  a  ■  Ta  —  Ke 


(28) 


By  integrating  (28)  from  0  to  H  or  z,  I<e  at  the  outlet  and  along 
the  z-axis  were  acquired.  Based  on  the  above  results,  the  distribu¬ 
tion  of  air  enthalpy,  air  humidity,  and  temperature  in  a  counter- 
current  adiabatic  dehumidifier  were  obtained.  Then,  the 
distribution  of  the  solution  parameters  can  be  calculated.  In  the 
paper,  the  analytic  solution  of  the  concurrent  flow  heat  and  mass 
transfer  process  were  also  given  out. 

Ren  et  al.  [52]  derived  a  new  analytical  solution  from  one¬ 
dimensional  differential  model.  By  introducing  some 
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Table  1 

Comparison  of  the  mathematical  models  for  adiabatic  dehumidifier. 


Classification 

Assumption 

Iteration 

Accuracy 

Applied  situation 

Finite  difference  model 

Least 

Extensive 

Best 

Component  design  and  operation  optimization 

e-NTU  model 

More 

Less 

Better 

Component  design  and  operation  optimization 

Simplified  models 

Most 

No 

Worst 

Annual  assessment 

Table  2 

Detail  information  of  the  mathematical  models  for  adiabatic  dehumidifier. 

Classification 

Model 

Flow  pattern 

Dimensionality 

Treatment  of  coupled,  heat  and  mass  transfer 

Finite  difference  model 

Factor  et  al.  [13] 

Counter 

One-dimensional 

Ackermann  correction 

Oberg  et  al.  [32] 

Counter 

One-dimensional 

Introduction  of  enthalpy 

Fumo  et  al.  [33] 

Counter 

One-dimensional 

Ackermann  correction 

Khan  et  al.  [34] 

Counter 

One-dimensional 

Ackermann  correction 

Liu  et  al.  [39] 

Cross 

Two-dimensional 

Introduction  of  enthalpy 

e-NTU  model 

Stevens  et  al.  [42] 

Counter 

One-dimensional 

Introduction  of  enthalpy 

Classification 

Model 

Flow  pattern 

Simplified  method 

Simplified  models 

Khan  et  al.  [34] 

Counter 

Correlations  based  on  simulation  results 

Liu  et  al.  [47] 

Counter  or  cross 

Correlations  based  on  experimental  results 

Gandhidasan  [48] 

Counter 

Introduction  of  dimensionless  parameters 

Chen  et  al.  [51] 

Counter  or  concurrent 

Introduction  of  parameters  similar  to  air  enthalpy 

Ren  et  al.  [52] 

Counter 

Introduction  of  dimensionless  parameters 

Babakhani  et  al.  [54] 

Counter 

Transfer  the  differential  equations  to  nonlinear  equations 

Liu  et  al.  [56] 

Cross 

Similar  method  of  cross-flow  heat  exchanger 

dimensionless  and  dimensional  groups,  the  conventional  equa¬ 
tions  of  one-dimensional  model  were  transferred  to  two  coupled 
ordinary  differential  equations,  whose  general  solution  are  as 
follows, 

AWM=C1ei'NTU‘  +  C2ei2NTU*  (29) 

A5=  -I<^C1e1'NTU‘  +  I<2C2el2NTUz  (30) 

The  model  is  just  suitable  for  the  case  where  the  solution  flow 
rate  and  concentration  change  slightly  as  it  assumed  that  the 
variation  of  the  equilibrium  humidity  ratio  of  solution  depended 
only  on  the  change  of  the  solution  temperature. 

Babakhani  and  Soleymani  developed  analytical  models  for  the 
counter  flow  adiabatic  regenerator  [53]  and  dehumidifier  [54],  The 
analytical  solution  was  deduced  on  the  basis  of  the  differential 
equations,  including  the  mass  balance,  air  humidity  and  tempera¬ 
ture  change  equations  listed  in  (1 )— (3),  and  liquid  desiccant 
temperature  and  concentration  change  derived  from  the  mass 
and  energy  conservation  equations.  To  achieve  the  simplification, 
two  main  assumptions  were  applied,  including  the  assumptions  of 
the  dilute  gas  phase  and  constant  equilibrium  humidity  ratio  on 
the  interface.  Then  the  above  equations  could  be  solved  and  the 
integrated  analytical  solution  were, 

W  =Wint+(Wi-Wint)exp(-aMNTUz)  (31) 

Tj,  =  Cj  +  C2exp(  -  0NTUz) - /  exp(  -  aMNTUz)  (32) 

(aMr  —  <xMQ 

T*  =  l? -  C2^xp(  -  emu^  +  jf11  n  eXP<  -  ^NTUZ)]  +  Ta 

RcLe  ( aM)z  —  aMO 

(33) 

lnX=  -Rm(Wi-Wint)exp(  —  aMNTUz)+C3  (34) 

Cs  =  Ga(Wj-Wint)exp(-aMNTUz)  +  C4  (35) 


With  the  above  solutions,  the  profiles  of  the  outlet  solution  and  air 
parameters  were  available. 

Based  on  the  above  model,  Babakhani  [55]  also  developed 
another  analytical  model  which  was  well  suited  to  the  high 
desiccant  flow  rate  conditions. 

Liu  et  al.  [56]  developed  an  analytical  solution  for  a  similar 
cross-flow  packed  bed  liquid  desiccant  air  dehumidifier,  whose 
numerical  model  had  been  reported  in  literature  [39],  In  the 
present  work,  it  was  regarded  that  the  desiccant  mass  flow  rate 
and  concentration  kept  constant  in  the  whole  process.  Thus  the 
Eqs.  (12)— (14)  were  got  rid  of,  and  only  Eqs.  (8)  and  (11)  were  left 
for  calculation,  which  were  rewritten  as 


dha 

H 

dhe 

0 

m* 

- —  ‘ 

-  = 

dZ 

L 

dX 

dha 

NTU 

(ha- 

-he) 

dZ 

=  ~ 

~1T 

It  was  found  the  above  control  equations  had  high  similarity 
with  those  of  the  cross-flow  heat  exchanger.  Thus  the  methods  in 
literature  [57,58]  were  referred  to.  As  the  solution  expressions 
were  a  litter  complicated,  here  they  will  not  be  presented. 

Recently,  Wang  et  al.  [59]  developed  a  simplified  yet  accurate 
model  for  real-time  performance  optimization.  Levenberg-Mar- 
quardt  method  was  used  to  identify  the  parameters.  The  proposed 
model  is  suggested  to  be  employed  in  real-time  performance 
monitoring,  control  and  optimization.  Park  and  Jeong  [60]  also 
developed  a  simplified  second-order  equation  model  as  a  function 
of  operation  parameters  to  study  their  impact  on  the  dehumidi¬ 
fication  effectiveness. 

3.4.  Summary 

To  sum  up,  there  are  three  models  to  simulate  the  adiabatic 
dehumidifiers.  The  general  comparison  of  them  is  presented  in 
Table  1.  As  a  matter  of  convenience,  the  information  of  some 
representative  models  is  summarized  in  detail  in  Table  2. 
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The  finite  difference  model  is  used  most  frequently  for  its 
accuracy.  However,  it  involves  complicated  iterative  process,  so  it 
is  only  suitable  for  the  component  design  and  operation  optimiza¬ 
tion.  The  common  assumptions  have  been  stated  in  Factor's  model. 
However,  some  other  simplifications  or  improvements  are  made  to 
satisfy  the  real  conditions.  It  can  also  be  observed  that  counter 
flow  configuration  is  the  most  commonly  used  flow  pattern,  and 
it  can  be  described  by  the  one-dimensional  model.  For  the 
dehumidifier  with  cross  flow  configuration,  the  problem  is  gen¬ 
erally  solved  by  the  two-dimensional  model. 

For  the  e-NTU  model,  two  additional  assumptions  are  included 
as  mentioned  above.  Thus,  the  model  is  less  accurate  than  the 
finite  difference  model.  However,  compared  with  the  finite  differ¬ 
ence  model,  the  calculations  are  dramatically  less  extensive.  Thus, 
the  e-NTU  model  has  great  potential  for  saving  time.  But  with  the 
development  of  the  computer  technology,  the  calculated  amount 
of  the  finite  difference  model  can  be  accepted  due  to  its  accuracy. 
Therefore,  much  fewer  researches  have  been  done  on  the  e-NTU 
model  in  the  later  studies. 

From  Table  2,  it  is  concluded  that  different  dimensionless 
parameters  have  been  introduced  in  the  deviation  processes  of 
the  simplified  models.  Meanwhile,  the  additive  assumptions  are 
needed  for  simplification.  Therefore,  they  are  only  applicable  for 
certain  operation  conditions,  and  different  models  can  be  chosen 
to  be  suited  to  the  real  situation.  The  biggest  advantage  of  these 
models  is  that  the  iteration  is  avoided.  As  a  result  of  their  high 
efficiency,  they  are  often  used  to  predict  the  annual  energy 
consumption  of  an  air-conditioning  system. 


4.  Models  for  internally  cooled  dehumidifier 

In  the  internally  cooled  dehumidifier,  some  cooling  fluid  is 
introduced  to  remove  the  heat  produced  by  the  vapor  condensa¬ 
tion,  as  shown  in  Fig.  7  [61  j.  Most  of  the  models  for  the  internally 
cooled  dehumidifier  are  developed  upon  the  finite  difference 
methods  used  for  the  adiabatic  dehumidifier.  The  difference  is 
the  additional  consider  of  the  heat  transfer  brought  by  the  cooling 
media.  In  addition,  because  of  the  relatively  low  solution  flow  rate 
in  the  internally  cooled  dehumidifier,  it  is  easier  for  the  solution  to 
form  a  thin  film  on  the  surface  of  the  padding  wall  or  plate.  Briefly, 
there  are  also  three  types  of  models  for  the  internally  cooled 
dehumidifier. 

Fluid  in  Solution  in  Air  out 

Gf,Tf  GS,TS,X  Ga ,  Ta+dTa 


Gf,T|+dTf  Gs+dGs  Ga,Ta 

Ts+dTs  ha,W 

XidX 

Fig.  7.  Heat  and  moisture  exchange  model  in  the  internally  cooled  dehumidifier 
[611. 


4.1.  Models  without  considering  liquid  film  thickness 

The  first  type  ignores  the  liquid  film  thickness.  Khan  and 
Martinez  [62]  developed  a  mathematical  model  to  predict  the 
performance  of  a  liquid  desiccant  absorber  integrating  indirect 
evaporative  cooling  to  achieve  an  almost  isothermal  operation. 
The  processed  air  and  the  solution  flowed  in  countercurrent 
direction  while  the  solution  and  water  flowed  in  parallel  direction. 


Primary  air 
stream 


Water  spray  Liquid  desiccant  spray 


Fig.  8.  Schematic  diagram  of  the  cross-flow  type  plate  heat  exchanger  [61]. 


a 

Centreline 


Desiccant 

solution 


Thin 

plate 


T 


dx 


Primary  air 


tV — 


b 

Centreline 


Fig.  9.  Schematic  diagram  of  the  control  volumes  considered:  (a)  primary  air- 
solution;  (b)  secondary  air-water  [63]. 
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Saman  and  Alizadeh  [63]  also  established  a  similar  model  for  a 
cross-flow  type  plate  heat  exchanger  (PHE)  serving  as  internally 
cooled  dehumidifier.  The  schematic  diagrams  of  the  PHE  are 
presented  in  Figs.  8  and  9.  It  is  obviously  that  it  has  the  same 
principle  with  that  in  literature  [62],  The  only  difference  lies  in  the 
configuration  difference.  Here,  Saman's  model  will  be  explained 
in  detail  as  it  is  more  representative.  In  the  paper,  the  PHE  was 
divided  into  a  certain  number  of  control  volumes  in  two  orthogo¬ 
nal  directions.  Several  assumptions  were  set,  including  no  heat 
transfer  with  the  environment,  negligible  temperature  gradient 
between  the  solution  film  and  water  film,  and  fully  cover  of 
solution  and  water  on  the  plate. 

Thus,  the  change  of  the  enthalpy  and  humidity  ratio  of  the 
primary  air  were, 


Desiccant  in  Air  out 


(Ts+dTs ,  X+dX)  (T,+dTa,W+dW) 


dha.p 

dx 


NTUs-Les 


H 


(he,s  ~  ha,p)  +  2  - '(We,s  —  Wp) 


(38) 


dWp 

dx 


NTUS 

H 


(Wes- 


Wp) 


(39) 


Also,  the  mass  and  energy  conservation  equations  in  the 
control  volume  were  given  as 


G' 


dWp 

dx 


(40) 


(41) 


Similarly,  the  change  of  the  enthalpy  and  humidity  ratio  of  the 
secondary  air,  and  mass  and  energy  conservation  equations  were 
expressed  as  follows 


dhq,r 

dy 


NTUf  ■  Lef 
H 


0h/~ ,r)"l~  yLe — (Wef  —  Wr) 


(42) 


dwr 

dx 


NTUf 

—R— 


■(We/-Wr) 


(43) 


dW 

dGr 

J  —  o 

dy 

dy 

dhar 

d(Gfhf) 

dy 

dy 

(44) 

(45) 


On  the  basis  of  the  above  governing  equations,  the  discretization 
equations  were  derived  for  each  control  volume.  The  calculation 
method  was  similar  to  the  flow  chart  of  the  finite  difference  model 
presented  in  Fig.  4,  while  it  was  more  complicated  as  a  result  of 
another  iteration  started  by  the  advance  assumption  of  the  cooling 
water  temperature.  In  the  paper,  it  is  noted  that  the  film  thickness  was 
mentioned,  but  the  simulation  did  not  take  it  into  consideration. 

Liu  et  al.  [21  ]  compared  the  performances  of  internally  cooled 
dehumidifiers  with  various  flow  patterns.  A  representative  heat 
and  mass  transfer  model  was  selected  for  detail  description,  as 
shown  in  Fig.  10.  Referred  to  literature  [34],  the  heat  and  humidity 
changes  of  the  air  were  almost  the  same  with  the  Eqs.  (6)  and  (8). 

The  energy  conservation  equation  was  expressed  as 


Cl 


dhq 

dx 


d(G'shs) 

dX 


+  CPjGf 


■  LdTf 
Hdy 


(46) 


As  the  mass  conservation  equation  has  been  explained  before, 
here  it  will  not  be  presented  for  limited  space.  And  the  heat 
transfer  between  the  desiccant  solution  and  the  cooling  water  was, 


djf 

dy 


NTUf 

L 


•( Ts-Tf ) 


(47) 


Combined  the  above  equations  with  the  inlet  conditions,  the 
distribution  of  the  parameters  of  the  air,  desiccant  solution,  and 
the  cooling  water  could  be  obtained.  The  model  was  applied  for 


Liquid  desiccant 

Humid  air 
- ^  Cooling  water 


Fig.  10.  Schematic  diagram  of  the  control  volume  for  a  counter-flow  internally 
cooled  dehumidifier  [21], 


analyzing  the  role  of  flow  patterns  in  depth,  which  was  seldom 
reported  in  previous  literatures.  In  the  later  study,  Zhang  et  al.  [64] 
employed  the  above  model  to  predict  the  performance  of  an 
internally-cooled  dehumidifier. 

Yin  et  al.  [25]  built  a  uniform  mathematical  model  for  an 
internally  cooled/heated  dehumidifier/regenerator  which  was 
made  up  of  a  plate  heat  exchanger.  It  is  important  to  point  out 
that  the  author  took  the  non-wetted  area  into  consideration 
by  introducing  the  wetness  coefficient.  Meanwhile,  in  a  control 
volume,  the  transfer  process  in  the  channel  width  was  considered 
to  be  symmetrical.  In  addition,  to  improve  the  accuracy  of  the 
model,  the  author  also  applied  the  correlation  of  the  mass  transfer 
coefficient  fitted  out  by  experiment  [24], 

Based  on  the  one-dimensional  differential  equations  for  the 
heat  and  mass  transfer  processes  with  parallel  or  counterflow 
configurations,  Ren  et  al.  [61]  developed  an  analytical  model  for 
internally  cooled  or  heated  liquid  desiccant-air  contact  units.  To 
increase  the  accuracy,  the  model  took  the  effects  of  solution  film 
heat  and  mass  transfer  resistances,  the  variations  of  solution  mass 
flow  rate,  non-unity  values  of  Lewis  factor  and  incomplete  surface 
wetting  conditions  into  consideration. 

Recently,  Qi  et  al.  [65,66  developed  a  simplified  model  to 
predict  the  performance  of  the  internally  cooled/heated  liquid 
desiccant  dehumidification  system.  Compared  to  previous  study, 
all  of  the  outlet  parameters  could  be  obtained  by  the  correlations 
quickly  and  accurately.  Thus,  it  is  very  useful  for  researching  the 
dynamic  operation  performance  of  the  internally  cooled/heated 
liquid  desiccant  dehumidification  system. 

Khan  and  Sulsona  [28]  developed  a  two-dimensional  and  steady- 
state  model  for  a  vapour  compression/liquid  desiccant  hybrid  cooling 
and  dehumidification  absorber  made  up  of  the  tube-fin  exchanger.  To 
simplify  the  model,  several  assumptions  were  made  reasonably.  The 
most  critical  one  was  to  consider  the  air  and  refrigerant  flow  in 
countercurrent  due  to  the  large  size  in  the  air  flow  direction.  In  this 
way,  the  complicated  problem  was  simplified  to  a  two-dimensional 
one.  The  governing  equations  are  very  similar  to  those  of  literature 
[21]  except  for  the  mass  conservation  equation.  There  were  three 
iterations:  one  for  the  real  refrigerant  exit  enthalpy,  one  for  the  correct 
local  coil  surface  temperature  and  the  final  one  for  the  actual  solution 
temperature. 

4.2.  Models  considering  uniform  liquid  film  thickness 

In  the  second  model,  the  film  is  considered  as  a  uniformly 
distributed  desiccant  film. 
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Park  et  al.  [26]  developed  a  three  dimensional  numerical  model 
for  simulating  the  coupled  heat  and  mass  transfer  in  a  cross-flow 
internally  cooled/heated  dehumidifier/regenerator.  The  schematic 
diagram  was  presented  in  Fig.  11.  Some  assumptions  were  used 
before  listing  out  the  governing  equations:  the  flow  was  consid¬ 
ered  as  laminar  and  steady:  the  physical  properties  for  both 
solution  and  air  were  constant;  species  thermo-diffusion  and 
diffusion-thermo  effects  were  negligible  and  thermodynamic 
equilibrium  existed  at  the  solution-air  interface.  Because  of  the 
relatively  small  absorption,  the  TEG  solution  film  thickness  was 
simplified  as  constant  and  the  film  mean  velocity  was  also 
unchanged.  Also,  the  velocity  gradient  in  the  liquid  solution  film 
at  the  interface  was  regarded  as  zero  and  the  flow  of  the  liquid 
solution  and  air  was  supposed  to  be  fully  developed  at  the  start 
place,  as  presented  in  Fig.  11. 

With  the  assumptions,  the  governing  equations  for  the  liquid 
solution  flow  were, 

cP"ll 

0  =  As“^2“+Ps&  (48) 


us^=Dt*l± 

dx  ts  dy2 


(49) 
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For  the  air  flow,  the  governing  equations  were  as  follows, 
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At  the  interface,  the  mass  and  energy  balances  also  existed, 
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The  above  equations  can  be  discretized  along  the  three  differ¬ 
ent  directions.  Being  combined  with  the  boundary  and  interfacial 
conditions,  the  equations  could  be  solved. 

Ali  et  al.  [67-69]  had  used  the  same  mathematical  model  to 
study  the  effect  of  the  flow  configuration,  the  inclination  angle,  the 
Reynolds  numbers,  various  inlet  parameters,  cu-ultrafine  particles 
volume  fraction,  and  thermal  dispersion  on  the  performance  of  the 
dehumidifier/regenerator. 

Mesquita  et  al.  [70]  compared  three  different  numerical  models 
for  parallel  plate  internally  cooled  liquid  desiccant  dehumidifier. 
The  second  one  introduced  a  constant  film  thickness.  It  was 
assumed  that  the  wall  was  isothermal  so  that  the  water  flow 
could  be  neglected.  On  the  basis  of  some  other  assumptions,  the 
dehumidification  process  could  be  described  by  a  two  dimensional 
model,  as  shown  in  Fig.  12.  Being  different  from  the  first  model 
which  does  not  consider  the  film  thickness,  this  model  took  the 
momentum  equations  into  consideration.  Firstly,  the  velocity 
profiles  of  the  air  and  liquid  desiccant  were  obtained  with  the 
momentum  equations,  presented  as  follows, 
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The  energy  and  species  equations  of  the  liquid  phase,  gas 
phase,  and  the  energy  and  species  balances  equations  at  the 
interface  are  almost  the  same  with  those  presented  in  literature 
[26],  Then  with  the  velocity  values  and  boundary  conditions,  all 
the  above  equations  could  be  solved  numerically  and  simulta¬ 
neously  with  the  software  package  Microsoft  EXCEL. 

Recently,  Dai  and  Zhang  [71]  employed  the  uniform  liquid  film 
thickness  to  evaluate  the  performance  of  a  cross  flow  liquid 
desiccant  air  dehumidifier  packed  with  honeycomb  papers.  The 
objective  of  the  paper  is  to  analyze  the  Nusselt  and  Sherwood 
numbers  in  the  channels  with  honeycomb  papers  as  the  packing 
materials. 


4.3.  Models  considering  variable  liquid  film  thickness 

The  final  model  introduces  a  variable  film  thickness.  Except  the 
constant  thickness  model,  Mesquita  et  al.  [70]  also  established 
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Fig.  12.  Schematic  diagram  of  control  volume  for  a  two-dimensional  model  [70]. 
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Table  3 

Detail  information  of  the  mathematical  models  for  internally  cooled  dehumidifier. 


Classification 

Model 

Flow  pattern  (air/desiccant) 

Flow  pattern  (desiccant/cooling  fluid) 

Cooling  fluid 

Regardless  of  film  thickness 

Khan  et  al.  [62] 

Counter 

Counter 

Water  and  air 

Saman  et  al.  [63] 

Counter 

Cross 

Water  and  air 

Liu  et  al.  [21] 

Six  different  configurations 

Water 

Yin  et  al.  [25] 

Cocurrent 

Cocurrent 

Water 

Khan  et  al.  [28] 

Cross 

Cross 

Ammonia 

Ren  et  al.  [61] 

Four  possible  flow  arrangements 

- 

Uniform  film  thickness 

Park  et  al.  [26] 

Cross 

Cross 

R22 

Ali  et  al.  [67-69] 

Cocurrent/counter/cross 

- 

- 

Mesquita  et  al.  [70] 

Counter 

Counter 

Water 

Dai  et  al.  [71] 

Cross 

- 

- 

Variable  film  thickness 

Mesquita  et  al.  [70] 

Counter 

Counter 

Water 

Hueffed  et  al.  [72] 

Cross 

- 

- 

Peng  et  al.  [73] 

Counter 

- 

- 

Diaz  [74] 

Cocurrent 

- 

- 

another  variable  thickness  model  for  internally  cooled  liquid- 
desiccant  dehumidifiers.  In  the  model,  the  film  thickness  in 
Eq.  (56)  varied,  thus  for  every  step  of  calculation,  the  film 
thickness  was  recalculated,  so  was  the  liquid  velocity  profile  in 
Eq.  (57).  However,  to  reduce  the  computational  time,  the  change  in 
the  airflow  velocity  profile  was  neglected  as  a  result  of  the  small 
film  thickness  changes.  Compared  the  results  of  variable  thickness 
model  with  the  constant  thickness  model,  it  was  found  that  the 
two  models  converged  to  the  same  results  at  higher  desiccant  flow 
rates.  However,  the  constant  thickness  model  underestimated  the 
dehumidification  for  low  desiccant  flow  rates. 

Hueffed  et  al.  [72]  presented  a  simplified  model  for  a  parallel- 
plate  dehumidifier,  with  both  adiabatic  and  isothermal  absorption. 
The  model  used  a  control  volume  approach  and  accounted  for 
the  film  thickness  variation  by  imposing  its  effect  on  the  heat 
and  mass  transfer  coefficients.  The  specific  equations  of  the  heat 
transfer  coefficient  in  terms  of  the  film  thickness  were  listed  as 


«c  = 


Nuka 


(60) 


dh  =  2(W-2S) 


(61) 


Then  the  mass  transfer  coefficient  was  obtained  from  the 
Chilton-Coulburn  analogy  as 


-2/3 


(62) 


In  each  control  volume,  various  parameters,  including  the  film 
thickness  S,  hydraulic  diameter  dh,  heat  transfer  coefficient  ac,  and 
mass  transfer  coefficient  aD  were  calculated  on  the  basis  of  the 
inlet  conditions.  In  this  way,  the  impact  of  the  film  thickness 
variation  on  the  heat  and  mass  transfer  process  were  introduced 
into  the  model. 

Recently,  Peng  and  Pan  [73]  investigated  the  transient  perfor¬ 
mance  of  the  liquid  desiccant  dehumidifier  with  a  one¬ 
dimensional  non-equilibrium  heat  and  mass  transfer  model. 
Unlike  the  previous  study,  the  local  volume  average  equations  of 
heat  and  mass  transfer  were  developed  in  the  work,  which  were 
solved  by  TriDiagonal-Matrix  Algorithm  (TDMA). 

In  the  later  year,  Diaz  [74]  also  developed  a  transient  two 
dimensional  model  for  a  parallel-flow  liquid-desiccant  dehumidi¬ 
fier.  The  geometrical  model  was  almost  the  same  with  that  of 
Mesquita  et  al.  [70],  The  difference  of  the  governing  equations  lies 
in  that  the  present  model  took  the  time  item  into  consideration. 
Some  non-dimensional  parameters  were  used  to  simplify  the 
calculation  process.  With  the  model,  the  variations  of  some  critical 


variables  over  time  were  illustrated  and  the  effects  of  oscillatory 
behavior  were  analyzed  in  depth. 

4.4.  Summary 

The  detail  information  of  the  mathematical  models  for  the 
internally  cooled  dehumidifier  is  listed  out  in  Table  3. 

The  models  for  internally  cooled  dehumidifier  without  con¬ 
sidering  the  film  thickness  are  very  similar  to  those  used  in  the 
adiabatic  dehumidifier.  However,  the  former  ones  are  more 
complicated  due  to  the  involvement  of  the  cooling  fluid.  These 
models  ignore  the  effect  of  the  velocity  field.  Thus,  the  results  of 
these  models  have  certain  discrepancy  with  the  reality,  as  velocity, 
mass  and  energy  have  strong  coupling  relationship. 

In  the  models  considering  uniform  film  thickness,  the  film 
thickness  and  velocity  are  usually  calculated  at  the  beginning 
and  then  keep  constant  through  the  whole  calculation.  The 
simulation  results  justified  that  the  constant  thickness  models 
under-predicted  the  dehumidification,  especially  for  low  desiccant 
mass  flow  rate. 

For  the  models  considering  variable  film  thickness,  all  of  the 
velocity,  mass  and  energy  equations  are  solved  together.  In  the 
process  of  each  iteration,  the  film  thickness  and  velocity  change,  so 
the  influence  of  the  flow  on  the  heat  and  mass  transfer  can  be 
demonstrated.  Thus,  the  model  is  most  accurate. 

5.  Conclusion 

Various  mathematical  models  have  been  proposed  to  assess  the 
performance  of  the  liquid  desiccant  dehumidifiers.  For  the  adia¬ 
batic  dehumidifier,  there  are  mainly  three  kinds  of  models:  finite 
difference  model,  e-NTU  model,  and  simplified  models.  The  finite 
difference  model  is  used  most  widely  for  its  accuracy.  However,  it 
involves  complicated  iterative  process  which  increases  the  com¬ 
puter  time,  so  it  is  only  suitable  for  the  design  of  the  component 
and  optimization  of  the  operating  conditions.  Like  the  finite 
difference  model,  the  e-NTU  model  also  requires  iteration,  but  it 
is  more  effective  and  less  accurate  correspondingly.  The  simplified 
models  require  more  assumptions,  so  they  are  only  suitable  in 
certain  operation  conditions.  But  due  to  their  high  efficiency,  they 
are  often  used  to  predict  the  annual  performance  of  the  system. 

For  the  internally  cooled  dehumidifier,  there  are  also  three 
kinds  of  models:  models  without  considering  liquid  film  thickness, 
models  considering  uniform  liquid  film  thickness,  and  models 
considering  variable  liquid  film  thickness.  The  first  model  ignores 
the  effect  of  the  velocity  field  totally.  Thus,  the  calculation  results 
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have  certain  discrepancy  with  the  reality.  The  second  one  consider 
the  effect  of  the  velocity  field,  but  the  assumption  of  the  steady 
state  and  fully  developed  field  cannot  describe  the  real  condition, 
especially  for  low  desiccant  mass  flow  rate.  The  final  model  is  most 
accurate  as  it  comprehensively  regards  the  influence  of  the 
changed  velocity  field,  but  the  calculation  becomes  relatively 
complicated. 

Though  a  large  amount  of  researches  have  been  conducted  on 
modeling  and  analyzing  the  liquid  desiccant  dehumidifier,  further 
efforts  are  still  needed: 

(1)  The  vast  majority  of  the  previous  models  assume  the  heat  and 
mass  transfer  process  to  be  steady  state  and  more  transient 
models  to  study  the  dynamic  performance  are  needed. 

(2)  More  three-dimensional  models  close  to  the  reality  are 
needed. 

(3)  Most  of  the  models  focus  on  the  outlet  parameters  and  more 
research  are  required  to  study  the  heat  and  mass  transfer 
process  in  the  dehumidifier  interior. 

(4)  It  needs  to  take  into  consideration  the  effect  of  variable 
physical  properties  which  are  taken  as  constant  in  almost  all 
existing  models. 
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